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We apply the optimal fluctuation method to the calculation of the optical absorption in disordered 
Oh' one-dimensional semiconductors below the fundamental optical gap. We find that a photon energy 

exists at which the shape of the optimal fluctuation undergoes a dramatic change, resulting in a 
different energy dependence of the absorption rate above and below this energy. In the limit when 
£ — | the interaction of an electron and a hole with disorder is stronger than their interaction with each 

other, we obtain an analytical expression for the optical conductivity. We show that to calculate 
the absorption rate, it is, in general, necessary to consider a manifold of optimal fluctuations, rather 
than just a single fluctuation. For an arbitrary ratio of the Coulomb interaction and disorder, the 
optimal fluctuation is found numerically. 
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INTRODUCTION 



I \ The interplay between disorder and inter-particle interactionsrfgsults in a number of remarkable phenomena, e.g., 
a singularity in the electron density of states at the Fermi energytlaBa and an enhanced localization length of pairs of 
interacting particles.Ej A familiar situation in which the simultaneous presence of Coulomb interactions and disorder 
plays an important role, is the process of optical absorption in disordered semiconductors below the fundamental 
' i' optical gap, i.e., below the band to band transition^ With the advent of modern optical materials, this problem 
is not only of interest in its -three-dimensional version, for whicK.it has received most attention, n but also in two 
dimensions (quantum wells )I3'B and one dimension (quantum wiresEro and semiconducting polymersa). Interestingly, 
in the absorption spectrum the relative strength of the Coulomb interaction between electron and hole and their 
interactions with the disorder, does not only depend on the disorder strength, but also on the photon energy: the lower 
' the photon energy, the larger the amplitude of a disorder fluctuation should be in order to create the corresponding 
Q\ , absorbing state below the fundamental optical gap. 

The effect of a relatively weak disorder on the Wannier exciton was considered in Refs. |l^ and [ll| Here, relatively 
weak refers to the situation where the exciton absorption peak is still visible as a separate peak below the band-to- 
band transition and is not smeared entirely by the disorder. Moreover, the restriction to weak disorder only applies 
for photon energies close to the exciton energy in the absence of disorder, E ex . In this case, the exciton localization 
length is much larger than the average electron-hole separation, and the exciton center-of-mass motion decouples 
from the relative motion of the electron and hole. The Wannier exciton then essentially behaves as a Frenkel exciton 
, in an effective disorder potential and the calculation of the optical absorption spectrum becomes a single-particle 
problem. This effective approach has been used to study numerically the_absorption and luminescence line shapes 
O , due to excitons in semi-conductor quantum wells with interface roughnessJij 

In Refs. [R] and [0], the low -energy tail of the exciton absorption p^&k was calculated analytically within the 
effective one-particle approach, by using the optimal fluctuation met hod. 11311-3 This method applies when the dominant 
contribution to the quantity of interest comes from disorder realizations close to a single large disorder fluctuation™ 
This method can be formulated as a saddle-point calculation of the functional integration oyer disorder realizations.^ 
It was used previously to calculate the optical absorption in disordered Peierls conductors^ and it is similar to the 
calculation of the Urbach tails resulting from the interaction of the electron-hole pair with the lattice, treated in the 
quasistatic approximation.!!!! 

In this paper we apply the optimal fluctuation method to the general case of arbitrary ratio of disorder strength 
and Coulomb interaction. In the strong disorder limit, when the exciton peak is destroyed, and for weak disorder, 
assuming sufficiently low photon energies, we obtain an analytical expression for the optical conductivity. We show 
that in these situations, the electron and hole in the optimal fluctuation are localized in two separate potential wells. 
We also show that in order to calculate the optical absorption rate, one has, in general, to consider a manifold of 
optimal fluctuations with different values of the electron-hole separation, rather than just one fluctuation. 

In the general case of arbitrary ratio of the Coulomb interaction and disorder strength the equation for the optimal 
fluctuation is solved numerically. We do this for an electron and hole described within the one-dimensional tight- 
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binding model. The equation for the optimal fluctuation m.this. model is similar to the equation for polaronic excitons 
and bipolarons obtained in the adiabatic approximation.l3EiE3 We show that, at some energy E c , lying in the region 
where the Coulomb interaction and disorder are of the same order, the shape of the optimal fluctuation undergoes a 
sudden change: while above E c it has the form of a single potential well, below E c it consists of a "dip" that localizes 
the hole and a "bump" that localizes the electron. This cross-over only takes place in one-dimensional systems. 

This paper is organized as follows. In Sec. II we introduce the continuum model that describes interacting electron- 
hole pairs in disordered semiconductors. To make the reader more familiar with the optimal fluctuation method, in 



Sec. Ill we briefly discuss its application to the calculation of the one-particle density of states. Then, in Sec. IV, 
we will discuss the shape of the optimal fluctuation for the electron-hole states. We argue that this shape crucially 
depends on the dimensionality of the system. We obtain and solve perturbatively the nonlinear nonlocal equation for 
the optimal fluctuation in the one-dimensional case. In Sec. [v| we obtain an expression for the tails of the optical 
absorption spectrum, which has a wider range of validity than the optimal fluctuation method. In Sec. [v| we present 
and discuss our results for the optimal fluctuation. Finally, we summarize and conclude in Sec. [VII] . Some technical 
details have been moved to appendices in order not to disturb the natural flow of the text. 



II. CONTINUUM MODEL AND OPTICAL ABSORPTION 



We consider direct gap semiconductors with Wannier excitons, in which case the low-energy electron-hole states can 
be described in the continuum approximation. In this approximation, the wave function W a (xi, X2) of the electron-hole 
pair with the energy E a , counted from the gap value A, satisfies the Schrodinger equation: 



2m/i 



Ai 



2m, 



■A 2 + tf(xi) - /3C/(x 2 ) + T/( Xl - x 2 ) 



(1) 



Here, xi and x 2 are the coordinates of, respectively, hole and electron, m^(m e ) is the effective hole(electron) mass, 
V(x) = — is the Coulomb interaction between the electron and hole, U(x) is the random potential due to impurities 
acting on the hole, while the disorder potential acting on the electron is —(3U(x). The dimensionless coefficient (3 
accounts for a different dependence of tho. energies of the bottom of the conduction band and the top of the valence 
band on the concentration of impurities lid We consider here the case of uncorrelated white noise disorder of strength 
A: 

(C/(x)C/(x')) =AS(x-x'). (2) 
The optical conductivity per unit volume for the electric field, polarized, e.g., in the x-direction, is given by 



a{u) = —Fihu - A), 



where the coefficient C, 
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C = 
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/ d 3 ru*(r)—v Q (r) 



(3) 



(4) 



is expressed through the periodic Bloch waves, uq(y) and i>o(r), describing, respectively, the electron and hole states 
with zero wave vector. In Eq.(Q) m is the electron mass in vacuum and the integration goes over one elementary unit 
cell Vb, in which the functions uq(t) and vq(t) are normalized to unity. 

The function F in Eq.(|^) is the part of the optical conductivity that has to be calculated within the continuum 
model: 



F(hu- 



A) = i ( \D a0 \ 2 S(A + E a - hoj) 



(5) 



In the last equation 



F>aO — 



I d d z**(x,x), 



(6) 



is the "continuum part" of the matrix element of the transition from the ground state to the excited state a (see 
Eq-(Q))- In Eq.(||) the brackets (. . .) denote the disorder average, V — L d is the total volume, L is the linear size, and 
d is the dimensionality of the system. 
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III. THE ONE-PARTICLE CASE 



In this section we briefly recall ho,™ thn optimal fluctuation method can be used to calculate the low-energy tail of 



the density states (per unit volume) 



13U14L1 



p(e) = ^(T l 6(e a [U\- e )Y (7) 

of a single particle moving in a random potential: 

/ ft 2 \ 

Hi> a = (-— A+£/(x)l ^ a =e a ^ a , a = 0,1,2,... (8) 
\ 2m J 

We are now interested in the density of states with a large negative energy. Such states can only be induced by 
large negative fluctuations of the disorder potential U (x) (in the absence of disorder the energy of all eigenstates is 
positive). The density of state is then, essentially, the probability to find such a fluctuation. When this probability is 
small, it suffices to keep in Eq. ([?]) only the contribution of the ground state (a — 0): 

p(e)K±(5(e Q [U\-e)), (9) 

because the probability to find a disorder fluctuation that induces an excited state with the same energy is even 
smaller. 

For the white noise potential Eq.(0) the disorder average can be performed by functional integration: 



p(e) = ~ JvUe-^5{e Q [U]-e) 



= 1 fouJ^e-Us+iWW-e)) (10) 
V I 2-niA 



S=\ Jd d xU 2 (x). (11) 



where the "action" S is given by 



The optimal fluctuation method is the saddle-point calculation of the functional integral in Eq. (|10|) , in which one 
assumes that the dominant contribution to this integral comes from the vicinity of one "optimal" disorder fluctuation 
C/(x), which has the highest weight among the disorder realizations that induce a state at the energy e. At the 
"saddle-point" the variation of 

S x = S + X(e [U]-e) (12) 

with respect to U(x) vanishes: 

The Feynman-Hellmann theorem allows us the carry out the variation of the energy with respect to the disorder as 
whence 

C/(x) = -A|^ (x)| 2 . (15) 

We must then insert this into the original Schrodinger equation (^) to find tpo- The resulting effective equation for ipo 
(which we can take to be real) is then the nonlinear Schrodinger equation 

h 2 

—Ai{jo + A^o + e^o = 0. (16) 
2m 
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It is convenient to introduce the dimensionless coordinate z = kx, where k is dchncd by e = — and the 

' J 2m 

dimensionless wave function d>(z): 



which satisfies 



A x <)> + 2<ff - 4> = 0. 



(17) 



(18) 



In dimensionless units the spatial extent of the wave function and the optimal fluctuation is of the order of 1. For 
d = 1, the solution of Eq.(18) can be easily found analytically 



1 



cosh z ' 



(19) 



while for d > 1, this equation has to be solved numerically. 
The solution of Eg. (|18|) is an extremum of the functional 



A[4>] 



d d z 



(V. 



a — b + c, 



(20) 



where we have used the notations 



a = jd d z(j) 2 , 

b = Jd d Z(/) i , 

{ c = fd d z(V^Y 



(21) 



The parameters of the optimal fluctuation can be expressed through a and b. Thus, the coefficient A in Eq. ( |l5| ) , found 
using the normalization condition Jd d x?pQ(x) = 1 and Eq.(|T7]), is expressed through a by 



A = a 



and the action S, defined by Eq.(|ll|), for the optimal fluctuation is given by 

h d \e\ 2 - d / 2 



s = 2 i-d/2 b: 



in 



d/2 



(22) 



One can obtain relations between the coefficients a, b, and c, using scaling arguments. In particular, if cf>(z) is a 
solution of Eq.(p"8|), then for the rescaled function 0a (z) = A</>(z), we have 



r/A 



A[(f>A 



0. 



A=l 



Since -A[</>a] = A 2 (a + c) — A 4 fo, we obtain 

a + c = 2b. 

Furthermore, the rescaling of coordinates, </>a(z) = 0(A _1 z) upon substitution into Eq.(|23|) gives 

(d-2)c + d(a - b) = 0. 

Combining Eqs.d24) and (Eq), we obtain 



(23) 

(24) 
(25) 



b = 



d 



a. 



We note that the same scaling arguments were used to prove that a scalar theory cannot have an inhomogeneous 
classical minimum in dimensions greater than 1 (the so-called Derrick's theoremcJ). Nonetheless, the inhomogeneous 



z 



FIG. 1: The shape of the wave function 4>{ z ) f° r d = 1,2, 3. 



solution 4>( z ) of Eq.(18) exists both for d — 2 and d = 3, as it gives an extremum of A (that is not bounded from 
below) rather than its minumum. 

The dependence of the function <fi on the radius z = |z| is shown in Fig. |l|. For d = 2, 3 we solved Eq.(|l^) numerically 
(see e.g., Refs. and while for d — 1 we used the analytical solution (|9|) for z > 0. The value of the coefficient 



a is 



2- 


for 


d = 


1. 


5.85... 


, for 


d = 


2. 


6.30... 


, for 


d = 


3. 



The result for the single-particle density of states at large negative energies, obtained by the optimal fluctuation 
method, has the form 

p(e)=K(e)e-*, (26) 

where the prefactor K(e) results from the Gaussian integrations over the small deviations of 8U(x) from the optimal 
fluctuation. The calculation of the prefactor also involves the integration over the locations of the optimal fluctuation, 
which cancels the volume V in the denominator of Eq.(Q). For d = 1 the prefactor K(e) is given byllj 

K{e) = ^. (27) 

An alternative derivation of this result, in which we use the correspondence between the optimal fluctuation and the 
instanton, describing the tunneling in the double-well problem, is given in Appendix A. 



IV. OPTIMAL FLUCTUATION FOR ELECTRON-HOLE STATES 



We now return to the problem of an interacting electron-hole pair in the presence of disorder. Similar to the previous 
section, we can obtain an equation for the wave function of the typical electron-hole state with large negative energy. 
However, before turning to a formal consideration (see Eqs. (B4T) and further), we first give a qualitative discussion 
of the properties of the optimal disorder fluctuation for the excitonic states. 

From Eq.((l|) we see, that the spatial extent of the disorder fluctuation that induces the typical single-particle state 
with the negative energy e, equals the spatial extent of this state, given by k _1 oc |e| -1 / 2 . Similarly, for the electron- 
hole state with negative energy E, the spatial extent of the optimal fluctuation, r oc |i?| -1 / 2 . While the magnitude of 
the disorder potential oc \E\, the magnitude of the Coulomb energy of the electron-hole pair oc oc IE] 1 / 2 . Thus, 
for large negative E, the Coulomb energy becomes smaller than the energy of the interaction of the electron and hole 
with disorder and the Coulomb interaction can be treated perturJaatively. This situation is very similar to that of the 
Coulomb gas, which becomes more ideal as its density increases.EJ 
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The shape of the typical electron-hole state crucially depends on the dimensionality of the system d. We first 
neglect the electron-hole interaction completely and consider a localized hole with the energy < and a localized 
electron with the energy e e < 0, such that Eh + £ e = E. Assuming that the electron and hole are localized far from 
each other, the action of the corresponding disorder fluctuation is given by the sum of the single-particle actions (see 
Eq.@): 



(I-. |2-d/2 I,- |2-d/2 
m h p^rrie 



(28) 



The factor /3 2 in the electron action is due to the fact that the strength of the disorder potential, acting on the 
electron, is 1 A. 

For d = 3 the minimum of the action is reached when only one particle is localized, i. e. at 

e h = E, e e = 0, for m h > m e (3 i / :i , 

e h = 0, e e = E, for m h < m e f3 4/3 . 

If we now include the Coulomb interaction between the electron and hole, we obtain the following picture of the 
typical electron-hole state with a large negative energy: One particle is localized by a disorder fluctuation, while the 
other particle forms a bound hydrogen-like state with the localized particle. If e.g. the hole is localized, then the 
total pair energy is 

E ~ e « + + (29) 

where = — h y Kh is the energy of the localized hole, the second term is the binding energy of the electron, 

e' 2 ) = — ^IfiS , and the third term is the energy of the repulsion of the electron from the disorder potential that 
localizes the hole. This picture is valid when \e^\ -C in which case the size of the hydrogen-like state, 

r = t^as (here as is the Bohr radius), is much larger than spatial extent of the disorder fluctuation localizing the 

hole kT 1 : 



WW 



In this limit the disorder potential acting on the electron can be approximated by 5— function and the corresponding 
energy, ~ 50 (jf^f") ^~Je\ — > ^ s parametrically smaller than \e^ 2%> \ and can be neglected. Then the absorption 



rate is: 



F(E) oc e = exp ■ 



bh 6 



2 1 / 2 ml /2 A 



E 



4 

m e e 



1/2- 



For d = 2, the situation is, essentially, the same as for d = 3: the minimization of the action (|2q ) gives that in the 
optimal fluctuation only one of the two particles is localized by disorder (the hole, for > f3^m e , or the electron, 
otherwise). For d = 1 the situation is quite different. In that case the minimum of the action (pq) is reached at 

^ = P 4 ^, (30) 

which means that in the one-dimensional electron-hole state with large negative energy both the electron and the hole 
are likely to be localized. In that case the optimal fluctuation has two parts (cf. Figs. [|(b) and ||(c) below): a part 
where the disorder potential is negative (to localize the hole) and a part with positive disorder potential (to localize 
the electron). The minimal value of the action is then given by: 



where M — rrih + (3 4 m e . 

There are two kinds of corrections to this action: (i) the correction due to the electron-hole interaction and (ii) the 
correction due to the interaction of the electron with the disorder fluctuation localizing the hole and vice versa. For 
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large negative E these corrections are relatively small. The regular way to find them is the perturbative solution of 
the equation for the typical wave function of the electron-hole state of a given energy. This equation can be obtained 
using the same considerations that led to Eq.(|l6|). The relation between the optimal fluctuation £/(x) and the wave 
function of the electron-hole state <3/(xi,X2) it induces, has the form 



C/(x) = -A J d d x' (* 2 (x,x') -/3* 2 (x',x)) , 



(32) 



where we used that \f r satisfies Eq.(|l|). Substituting this optimal fluctuation back into the Schrodinger equation (|l|) 
leads to a nonlinear, nonlocal equation for \E': 



2m, 



-Ai 



2m P 



A 2 + V(x 2 - xi) - E U( Xl , x 2 ) = A*(xi, x 2 ) 



< jd d x' (* 2 (x!,x') - /3* 2 (x', Xl ) - /3* 2 (x 2 ,x') + / 3 2 * 2 (x',x 2 )) 



(33) 



ion for the excitonic polaron and bipolaron wave functions, obtained 
>0 ' It is clearly impossible to solve this equation analytically. In the 



This equation is formally equivalent to th- 
in the adiabatic treatment of the lattice 
remainder of this section we obtain the action of the optimal fluctuation for d = 1 using a perturbative solution of 
this equation, while in Sec. VI we give the results of a numerical solution, also for d = I . 
The first term in the perturbative expansion of the wave function, 

* = * + *i + • ■ • , 

is the product of the wave functions of a noninteracting electron-hole pair separated by some distance: 



^o(xi,X 2 ) = 1ph(xi - Xh)lpe(x2 " X e ) = — £ " (f){K h (xi - Xh))<f>(K e (xi - X e )), 



(34) 



where Xh{x e ) are the average hole(electron) positions, the single-particle wave function <fi is given by Eq.(U9) and the 

wave vectors Kh and K e are defined by = — an d £ e — — -^r ■ 

In the zeroth order of the expansion the disorcter potential in the optimal fluctuation is obtained by substituting 
Eq.(pl) into Eq.(||): 



U (x) = U h (x - x h ) + U e (x - x E ) = -\ipl(x - x h ) + j3X^ 2 e (x - x e ), 



(35) 



where the first term is the negative potential (a dip) that localizes the hole and the second term is the positive 
potential (a bump) that localizes the electron. 

Inserting Eqs.(|3J) and ( |35| ) into Eq.(|3"3|) and taking into account that £/, + e e — E, we obtain that the cancellation 

of the zeroth-order terms in that equation requires A = 2h Kh and A/3 2 = 2h K ° , from which we find 



m h zp 
At ' 



M 



E, 



(36) 



2\E\ 
M 



A = 2hy 

which agrees with Eq. (^p|) obtained above. Also, the calculation of the action in the zeroth-order approximation, 



dx[Ul(x) + U?(x)] 



Sn = 



gives Eq.(|3_y). 

Though the calculation of the first-order correction to the wave function is, in general, difficult, the first-order 
correction to the action can be expressed through the unperturbed wave functions of the electron and hole: 



Si (r) = A / dxtp h (x) 



dx'ip 2 (x')V(x -x' + r) + (3\tp 2 {x - r) 



(37) 
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where r = x e — Xh is the average electron-hole separation. The somewhat lengthy derivation of this result is given in 
Appendix ||. The first term in the square brackets is the correction due to the electron-hole interaction, while the 
second term describes the interaction of the hole with the bump, localizing the electron and vice versa (see Eq.([35|)). 
The optimal electron-hole distance r* is found by minimizing S\ with respect to r: 



dSi 
dr 



0. (38) 



From the form of Eq.(p7[) it is clear that the latter condition is just the balance of average forces acting on the 
hole: the attraction from the electron and repulsion from the disorder potential localizing the electron. Thus the 
typical electron-hole state of large negative energy can be considered as a kind of "molecule" , in which the disorder 
fluctuations binding the electron and the hole play the role of "nuclei" with, respectively, positive and negative charge. 



The result Eq.(37) can be cast into a more transparent form using the fact that for the Coulomb interaction between 
the electron and hole to be small compared to their interaction with the disorder fluctuations, the optimal electron- 
hole separation r* should be large compared to the spatial extent of the disorder fluctuations: eT KhT * , e~ K " r * -C 1. 
We furthermore assume, that mh > m e (3 4 and that e -( K h-«e)»'» <g; i ; j. e-) the hole is sufficiently more localized than 
the electron, which holds unless nih is very close to m e (3 4 . Then Eq.pTn becomes 



A er 

For the optimal distance r», at which Si has its minimum, we obtain 



S ' (r) ~ -+ml(r). (39) 



-Vs** (40) 

where A is given by Eq.([j6|). The dimensionless parameter on the right-hand side is, essentially, the square root of the 
ratio of the Coulomb energy to the total energy E, which, by assumption, is small (and thus K e r* is logarithmically 
large). 

Furthermore, the second term in Eq.(|3^) (due to repulsion of the hole from the electron optimal fluctuation) is 
small compared to the first term (due to the electron- hole interaction): 

ml{r*) 1 



e 2 /(er*) 2K e r» 

Therefore, 



S ^-x—^^ — 

er* dE er* 



where in the last step we used 



as follows from Eqs.(|3l|) and (|3^). We thus see that, to the lowest order, the effect of the correction to the action Si 
is to replace S (E) by S (E + -^-): 



F(£)ocexp<f 4 ^ 



3AM 1 / 2 



E+ — 



where the electron-hole Coulomb shift depends on E, as according to Eq. (fio|) , r* cx |i?| -1 / 2 In \E\. 

To end this section about the shape of the optimal fluctuation we note that the numerical solution of Eq. (^) for 
sufficiently small \E\ gives r* = 0, which corresponds to the localization of the electron- hole pair by a symmetric 
single-well disorder potential. In that case the repulsion of the electron from the disorder potential that localizes the 
hole is compensated by the Coulomb interaction between the electron and hole. The transition from the symmetric 
to the asymmetric shape of the optimal fluctuation is studied in detail Sec. refnumres, where the results of the 
numerical solution of Eq.(^) are discussed. We shall show that our analytical approach, in fact, gives a rather 
accurate description of that transition (see Fig. ^|). 
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V. OPTICAL ABSORPTION TAIL IN ONE DIMENSION 



In the previous section we have shown, that for d = 1 the electron and hole in the typical state with a large negative 
energy are localized relatively far from each other. This property allowed us to calculate approximately the weight 
of the optimal disorder fluctuation for the electron-hole pair. It also simplifies the calculation of the pre-exponential 
factor in the expression for the absorption rate, which results from the integration over small deviations from the 
optimal disorder fluctuation. In this section we obtain an analytical expression for the tail of the optical absorption 
spectrum, which actually has a wider range of validity than the standard optimal fluctuation method. 

As we have shown in the previous section, the dominant contribution to the optical absorption at large negative 
E comes from the disorder realizations that are close to the sum of the single-particle optimal fluctuations, Uq(x) = 
Uh(x — xh) + U e (x — x e ) (see Eq.(^5|)), where Uh(x — x e ) localizes the hole with energy en near x — Xh and U e (x — x e ) 
localizes the electron with the energy e e near x — x e ■ In the previous section we have calculated the optimal electron- 
hole separation, for which the weight of the disorder fluctuation reaches its maximum. In this section we shall treat 
Xh and x e , as well as the single-particle energies Sh and s e , as the collective variables of the functional integration 
over disorder. Performing the Gaussian integration over all the small deviations SU(x) = U(x) — Uq(x) in Eq.(||) that 
are orthogonal to the deviations corresponding to the four collective coordinates, we obtain 

F(E) * iJde h de e dx h dx e K h (s h )K e{ s e )D^x e -x h ) 

x e~i {Sh{£h}+s " i£ ' }+ss) S{e h + s e + 6e - E). (42) 



Here, L is the chain length, D is the transition matrix element (see Eq.(|6|)), Sh(£h) and S e (e e ) are the single-particle 
actions, given by Eq.(E8|), and SS = SS(r), where r = x e — Xh, is the correction to the action of the electron-hole pair 
due to the overlap between the electron and hole optimal fluctuation: 

8S(r) = JdxU h (x - x h )U e (x - x e ). 

Furthermore, Se — 5e(r) is the energy correction due to the electron-hole interaction and the interaction of the 
electron with hole optimal fluctuation and vice versa. The first-order correction, calculated using the unperturbed 
electron- hole wave function Eq.(R3), is 



Se(r) = J dx [ij%(x — Xh)U e (x — x e ) + ipl(x — x e )Uh(x — Xh)] 

+ dxdx'tpl(x - x h )V(x - x')il? e {x' - x e ), (43) 



where r = x e — Xh is the electron-hole separation. In the same approximation the transition matrix element is given 
by 



D(r) — Jdxip h {x - x h )ipe{x - Xe)- 

Finally, in Eq.([f2"|) we used that for small overlap between the electron and hole optimal fluctuations the prefac- 
tor, resulting from the Gaussian integration over 8U(x), is the product of prefactors for isolated hole and electron, 
K h (£ h )K e (£ e ))El (We note, that in K e {s e ) the disorder strength A has to be substituted by f3 2 A.) 
The integration over the center-of-mass coordinate 

m h x h + m e x e 



R 



run + m e 



which is a zero mode, gives the chain length L. Due to the <5-function in Eq.([42|) e e = E — Se — eh and the remaining 
integration over the hole energy eh can be performed in the saddle-point approximation. The condition of the minimum 
of Sh(eh) + S e {e e ) gives Eq.(|30|), which determines the saddle-point values of eh and e e . Then Eq.(|4^) can be written 
as follows: 



B 

where 



F(E) ~ I drK h (s h )K e (e e )D 2 (r)e-^ E -^+ ss \ (44) 
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FIG. 2: The energy dependence of the logarithm of the absorption rate. 



_ m c m h 



and S (E) is given by Eq.(|l|). 

Expanding S (E-Se(r)) w S (£) - f| fe(r) and using Eq. (^lj) and the relation between the single-particle optimal 
fluctuations and the wave functions: Uh(x — x^) = —Xi/j^x — x^) and U e (x — x e ) = Xf3ip 2 (x — x e ), we find that the 
r-dependent part of the action coincides with the first-order correction to the action of the optimal fluctuation Eq. (^7|) : 



5S- 



Thus, we can write Eq.(B3) in the form 



F{E) 



yj^- J drK h {s h )K e {e e )D\r)e-^ s ^ s ^\ 



(45) 



The typical energy dependence of the absorption rate Eq. ([l5|) is shown in Fig.|[ where we plot log F as a function 
of (E I Eq) 3 ^ 2 , with Eq < being the exciton binding energy. One can see, that the energy-dependence of logF very 
quickly becomes linear for energies below Eq. The linear dependence reflects the relative weakness of the Coulomb 
interaction between the electron and hole (see Eq.(31)), which was the main assumption of this analytical calculation. 

In the case when the integral over the electron-hole separation r comes from a small vicinity of the saddle-point 
r*, determined by Eq.(38), we find that the saddle-point action coincides with the action for the optimal fluctuation, 
So + Si(r«), obtained in the Sec. [V. The result of the saddle-point integration over r is 



F(E) 



2ttA 



D 



iir- 



Ph(£h)p e (£e)D 2 (r*)e" 



where 



d 2 S 
dr 2 



2Ae 2 



(K e r» - 1) 



Equation (45) also applies when 
realizations t 



and eh and e e are the unperturbed single-particle energies, given by Eq.(^) (we have used that the factors Kh(£h) 
and K e (e e ) are slow functions of the energies). 

instead of one optimal fluctuation, one finds an entire manifold of disorder 
lat significantly contribute to the optical absorption. In particular, for a non-interacting electron-hole 
pair the minimum of the action Si is reached at the largest possible electron-hole separation, since for V(xi — X2) = 
nothing can counteract the "repulsion" between the two disorder fluctuations, described by the second term in Eq.(|37|). 
However, these electron-hole states clearly do not contribute to the optical absorption as the transition matrix element 
D vanishes at infinite r. Since the repulsion decays exponentially with the electron-hole separation, the action S\ is a 



very weak function of r, as soon as the latter exceeds the spatial extent of the electron and hole states, 



and k„ 



The contribution of the electron-hole pairs with large r is then suppressed not by the weight of such fluctuations, but 
by the smallness of the transition matrix clement D, which decays exponentially with r. Thus, when the electron-hole 
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interaction is absent or relatively weak, to calculate the absorption rate, one has to sum the contributions of many 
disorder fluctuations with different electron- hole separations, which can be accomplished using Eq.(|45|). 

The integration in Eq.(|4|) for non- interacting electron and hole becomes particularly simple in the case when 
the hole is localized stronger than the electron, or more precisely, when e ~( K ' l ~ fie ) r <C 1 for relevant electron-hole 
separations r, which we already have used in the Sec. IV. Then the action Si for the non-interacting case is 



Si(r) = X 2 f3 Jdxilj 2 h (x)iJ 2 e {x + r) » \ 2 {3ip 2 e (r) 
and the transition matrix element is given by 

D(r) 

and 



/ 2~Kh 



drD{r) 2 e 



Sl(r-) 



7-oo cosh 2 z 



dz < 

e cosh 2 



TT 2 0MA 

32 M |£| 2 ' 



where we assumed that 



, \ 2 I3k £ 3m e f3 3 So 
C = ~ . = — 7T— -r > 1. 



2A 



M A 



(46) 



Since for the applicability of the optimal fluctuation method, anyhow, the action So has to be much larger than A, 
Eq.(E3) holds, unless 3m jf is very small. 



Thus, finally, for non-interacting electron-hole pairs the absorption rate is given by 



1 


7TjU 


l\E\ 


4 I 


~ 


2Ah\ 


1 2M 


exp < 



4h 2\E[ 



3A\ M 



VI. NUMERICAL RESULTS 

The analytical results of sections IV and for the optimal fluctuation and the optical absorption spectrum were 
obtained by a perturbative treatment of the Coulomb interaction between the electron and hole. When this interaction 
is of the same order as the magnitude of the disorder potential the optimal fluctuation has to be solved numerically. 
In this section we present our numerical results for the optimal fluctuation in one dimension. 

Instead of solving Eq.(|33|) directly, we perform the numerical calculations for a tight-binding model defined on a 
lattice. The discrete version of the Schrodinger equation ([!]) reads: 

- ti{lpn-l,m + 1pn+l,m) - t2(lpn,m+l + V'n.m-l) + V nm lp nm + (U n - (3U m )lp nm = Slpnm, (47) 

where t\ and t% are, respectively, the hole and electron hopping amplitudes, the indices n, m = 1, 2, . . . , L denote the 
sites of the one-dimensional lattice with the lattice constant a and periodic boundary conditions, U n is the disorder 
potential: 

{U n ) = 0, {U n U m ) = — S nm , 
a 

and V nm is the regularized Coulomb interaction 

V n-m = -90[-\ r~T~T + i T~T )• (48) 

\\n-m\+o n ,m Jy - \n - m\ + o\ n - m \,N J 

2 

Here, go — — and the second term in the brackets is added to satisfy the periodic boundary conditions. 

The discrete equation ( |47| ) reduces to the Schrodinger equation (ID) in the continuum limit, when all relevant electron 
and hole states have small wave vectors, k±a, <C 1. In that case, the dispersion of the free hole dispersion is 

£h{k) = —2t\ coska w —2t\ + ti{ka) 2 , 
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so that rrih — 2 t 1 a i anci ) similarly, for the electron we have m e — 2t " a ^ . In this calculation we put a = t\ = t% = 1. 
Then the continuum limit is reached for small values of the coupling constant go and the energy E, counted from the 
bottom of the band: 



E = £ + 2(tx+t 2 )=£ + 4<l. 
The discrete analog of the action ( |l2| ) is 

SxW,U] = 5^Dj + A»^-£), (49) 



where 



r , _ — 2Xi nm [V , nm(*lV'n+l,m + *2 ) + V nm ipl m + (U n - i3U m )ipl m ] 
&{U,ty\- =; • (DUJ 



The denominator in the last equation takes care of the wave function normalization. It is readily seen that varying 5a 
with respect to the wave function ip nm (that can be chosen real) yields the discrete Schrodinger equation (47), while 
the minimization of S\ with respect to the disorder potential U n gives a relation between the optimal fluctuation and 

U n = -TrY,Wn™-M 2 mn), (51) 



where we have used the notation 



AA = 5>L- (52) 



We find the optimal fluctuation U n and the corresponding wave function ip nm by minimizing the functional 

2 \ 



nm nm 



n \ m 




(53) 

which is the two-particle analog of Eq.(|2C|). While S\ depends on both ij) nm and U n , A\ is a functional of the wave 
function only. One can easily check, that the condition ^ x = is equivalent to Eq.(|47j) with the disorder potential 



U n given by Eq.(51). The minimization of A\ with respect to ip nm was carried out numerically, using the steepest 
descent algorithm for L = 50. The energy of the electron-hole pair and other quantities of interest, e.g., the optimal 
fluctuation U n , its weight, and the corresponding electron-hole wave function, are first obtained as functions of the 
Lagrangian multiplier A. Then we eliminate A by replotting these quantities as functions of the energy E. 

We first consider the energy dependence of the weight of the optimal fluctuation, which, essentially, determines the 
energy dependence of the optical absorption rate. This weight is given by e - ^, where S = | J2 n ^ s the action of 
the optimal fluctuation [cf. Eq.(|Tl|)]. Motivated by the dominant \E\ 3 / 2 behavior of S [cf. Eq.([5l])], we plot in Fig. || 
S 2 / 3 as a function of the energy E for go = 0.2 and (3 = 0.5. The open circles are obtained by the numerical procedure 
described above, while the solid line is the result of our approximate analytical calculation of the action: S = Sq+Si, 
where So and Si are given by Eqs. ( |3l| ) and (|37|), respectively. Clearly, apart from a small energy interval near the 
exciton binding energy in the absence of disorder, Eq w —0.08, the energy dependence of S 2 ^ 3 is indeed close to linear 
and the agreement between our numerical and analytical results is good. 

For energies close to Eq, the assumption that disorder dominates the Coulomb interaction, used in our analytical 
approach, breaks down, which explains the deviations of the numerical data from the S oc |£?| 3 / 2 law and from the 
analytical curve. On the other hand, the deviations found at relatively large energies \E\ ~ 1 are due to the break 
down of the continuum approximation, resulting from the fact that the hole becomes localized on a single lattice site. 

These changes in the energy-dependence of S reflect changes in the shape of the optimal fluctuation and the 
corresponding electron-hole wave function. In Fig. ^ we plot U n and ip nm , calculated numerically, for three different 
values of the energy: —0.14, —0.30, and —0.94. The first value E = —0.14 is rather close to the exciton binding 
energy Eq. In that case, the optimal fluctuation, shown in Fig. ^a, is rather shallow and it is symmetric around 
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FIG. 3: The action S of the optimal fluctuation to the power 2/3 plotted as a function of the electron- hole energy E for 
go = 0.2, f3 = 0.5 and the chain length L — 50. The open circles are results of the numerical calculation and the solid curve 
was obtained analytically (see explanations in the text). 



xq, where xq is the position of the minimum of this fluctuation. This symmetry implies that no separation exists 
between the average electron and hole positions. Such an optimal fluctuation was discussed in Ref. |ll| for the situation 
where the Coulomb interaction dominates the disorder and the spatial extent of the optimal fluctuation l opt is much 
larger than the exciton radius r exi leading to decoupling of the center-of-mass and the relative motion. That limit 
is rather difficult to simulate numerically within our discrete model, as the requirement to maintain the validity of 
the continuum approximation then leads to 1 <C r ex <C l opt <C L, thus forcing us to consider a very large lattice size 
L. From Fig. |^ one observes that for E = —0.14, l op t is comparable to the exciton radius r ex ~ 5. Still, one can see 
from Figs. 0d and g, that the electron-hole wave function, apart from the derealization along the electron coordinate 
X2, also shows a strong derealization along the line x\ = X2, which corresponds to the center-of-mass motion of the 
exciton. For the second value of the energy, E = —0.3, the optimal fluctuation has the asymmetric "dip-bump" shape, 
which corresponds to the localization of the electron and hole on different sites of the chain (see Fig. [|b). Finally, at 
E = —0.94, the hole is practically localized on one chain site, in which case the discrete model should not be used to 
simulate the continuum one (see Fig. ^c). From Figs. f| e, f, h, and i, one can see that for E = —0.3 and E = —0.94 
the electron-hole wave function tb nm is mostly delocalized along the electron coordinate x-i and it has no derealization 
along the center-of-mass direction, x\ =x%. 

The cross-over from the symmetric to the asymmetric shape of the optimal fluctuation can be most clearly seen 
from the energy dependence of the average electron-hole separation, defined by 

nm 

The results of the numerical calculation of r as a function of the energy E are shown in Fig. || by circles. In that 



figure we also plot the optimal electron-hole distance r*, obtained by our approximate analytical approach of Sec. IV 



[Eq.(38)]. The minimization of the correction to the action S\ was performed numerically, as the approximation that 
was used to obtain Eq.(^0|) is too crude to describe the changes in the shape of the optimal fluctuation. Figure || 
shows a very fast transition from the optimal fluctuation with zero average electron-hole separation to one with finite 
separation. One can also see that the approximate analytical approach provides a good qualitative description of this 
cross-over (in particular, the onset of the cross-over and the shape of the r(E) curve), but it gives a somewhat larger 
value of the electron- hole separation at low energies. 

To clarify the nature of this cross-over, we plot in Fig. ||(a-c) the coordinate distributions of the electron (solid line) 
and hole (dotted line), 



2 



P h = V ib 

p e = y ib 2 

for g — 0.2, j3 — 0.5, and three different values of energy E. In addition, Figs. ||(d-f) show the corresponding effective 
Hartree potential acting on the electron 



14 




(c) 


1 


I 



25 
X 



50 




FIG. 4: The shape of the numerically obtained optimal fluctuation for E = —0.13 (a), E = —0.3 (b), and E — —0.94 (c). 
Panels (d - f) show the corresponding electron-hole wave function tp„m and panels (g - i) show the contour plot of the wave 
function. All plots correspond to go = 0.2, j3 = 0.5, and L = 50. The coordinates x\ = na and X2 = ma, where n,m — 1, . . . , L 
describe, respectively, the hole and electron positions in the chain in units of a = 1. Note, that as E decreases, the shape of 
the optimal fluctuation undergoes a transition from a single-well to a "dip-bump" structure. 



At E — —0.18, just above the cross-over, the electron wave function has two peaks, in accordance with the double-well 
structure of the effective Hartree potential (see Figs. || a and d). The peak separating the two potential wells is the 
disorder potential that localizes the hole and repels the electron. As the energy E decreases, the height of the peak 
grows, which suppresses the electron tunneling between the two wells. For large separations between the wells, the 
weight of the symmetric optimal fluctuation, in which the electron is delocalized over the two wells, is lower than the 
weight of the single- well optimal fluctuation with the same electron energy. For E = —0.20 in the cross-over region 
the two wells become unequal [see Figs, and |^e], and at E = —0.50, well below the cross-over region, the electron 
is predominantly located in a single well [see Figs. ||c and^f]. 



VII. CONCLUSIONS 



In this paper we studied theoretically the photoexcitation of electron-hole pairs in disordered one-dimensional 
semiconductors. Using the optimal fluctuation method, we calculated the low-energy tail of the optimal absorption 
spectrum in these systems. We were in particular interested in the effects of the Coulomb interaction between the 
electron and hole on the energy dependence of the absorption spectrum. We want to point out, however, that the 
calculation of the absorption rate is a nontrivial problem even for non-interacting electron and hole. In particular, it 
cannot be reduced to a single-particle calculation, since, on the one hand, the electron and hole move in a common 
disorder potential, and, on the other hand, the effect of disorder on the electron is different from the effect of disorder 
on the hole. 

We found that, as the photon energy decreases, the shape of the optimal fluctuation undergoes a cross-over. Close 
to Eq (the exciton binding energy in the absence of disorder), the Coulomb energy dominates disorder and the optimal 
fluctuation has a symmetric shape. It reflects the fact that at those energies the exciton is not entirely destroyed 
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FIG. 5: Numerical (circles) and analytical (solid line) energy-dependence of the average electron-hole separation calculated 
for g = 0.2, (3 = 0.5, and L = 50. 
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FIG. 6: The left-hand side of the figure shows the electron (solid line) and the hole (points) distributions in the typical 
electron-hole state for E = -0.18 (a), E = -0.20 (b), and E = -0.50 (c). On the right-hand side (panels (d) - (f)) the 
corresponding effective potentials acting on the electron are plotted. All plots correspond to go — 0.2, = 0.5, and L = 50. 

by disorder. That limit was considered previously in Ref. [ll]. In the opposite limit, when the disorder dominates 
the Coulomb interaction, the optimal fluctuation has two parts: a "dip" that localizes the positively charged hole 
and a "bump" that localizes the negatively charged electron. This cross-over to an asymmetric optimal fluctuation 
is characteristic for motion in one dimension - it does not occur in two and three dimensions. This cross-over can, in 
principle, be observed experimentally, as the absorption rate has a different energy-dependence above and below the 
cross-over energy (see inset in Fig. |3j) . 

At photon energies well below the exciton energy, when the Coulomb interaction can be treated as a perturbation, 
we obtained an analytical expression for the optical conductivity of a disordered one-dimensional semiconductor. We 
showed that in some cases one has to go beyond the standard optimal fluctuation method and consider an entire 
manifold of optimal fluctuations, corresponding to different electron-hole separations. We also performed numerical 
calculations of the optimal fluctuation for a discrete model that allowed us to study the whole region of photon 
energies, in which the optimal fluctuation method is applicable (i.e., for E < Eq, for which the action S is much 
larger than A). In the continuum limit we found a good agreement between our numerical and analytical results. 
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APPENDIX A: SADDLE-POINT INTEGRATION 



The integration over disorder realizations close to the optimal one for a particle in a one-dimensional white-noise 
potential was performed by a number of different methods. Ej Here we perform this integration by reducing it to a 
path-integration over the trajectories of a particle moving in the double-well potential. A more detailed discussion of 
this method for a similar problem is contained in Ref.E3. 

To find the average density of states, we have to perform the functional integration 



L 



VU 



+ L/2 
L /2 



dxU* 



(Al) 



In this appendix we put = 1, so that Eq. (0) becomes 



E a 1p a (x), 



(A2) 



and e = —n 2 . We consider Eq.(A2) on the interval x S [— i, + -§■] . 

The correspondence between the disorder average and the quantum-mechanical average is established through the 
functional substitution U(x) — > y{x): 



U(x) 



dy 

dx 



y 2 +e. 



(A3) 



Then the action for the disorder realization U(x) can written in the form 



s- 1 - 



r +L/2 




/ dxU 2 (x) 


l-L/2 




,+L/2 

/ dx 


'(dy 


l-L/2 


\dx 



K 2 y ■ 



y_ 

3 



+L/2 
-L/2 



(A4) 



The integral in the second line of Eq.( A4) can be considered as the Euclidean action of a particle with the coordinate y 
moving in the double-well potential W(y) — ^ (y 2 — k 2 ) . Thus, the functional integration over disorder realizations 
U(x) can be represented as the path-integral over trajectori es o f the particle y{x) in the imaginary time x. 

The regularized form of the functional substitution Eq.(A3), obtained by discretizing the imaginary time x n 
— 4 + an, where a = -4 and n — 0, 1, . 



,N, is 

(y n -y n -i) 



a 2 

Vn-l 



e-U n 



ay n-i 



Here y n = y(x n ), U n — U(x n ), and terms of order a 2 are neglected. The last equation is a one-to-one correspondence 
between U (x) and y{x) if the value of y at the left boundary of the interval of x is kept fixed: 

2/o = y(~L/2) = y L , 

where yL is an arbitrary number. The Jacobian of this functional substitution is given by 



J = 



VU 



exp 



+L/2 



L/2 



dxy(x) 



Furthermore, if, for a given U(x), e — e a is the energy eigenvalue, then y(x), defined by Eq.(A3), is related to the 
eigenfunction ip a (x) by 



2/0) 



In that case y(x) also satisfies a condition at the right boundary +j± 



lpa(x) ' 

idar 

UN = y(+L/2) = y R . 
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Therefore, the density of the single-particle states can be written in the form 

dy(L/2) 



ds 



S(y(L/2)-y R ). 



where the derivative ^ at x = +L/2 is calculated at fixed y(—L/2) and U(x). Differentiating Eq. (|A3|) with respect 
to e, we obtain 



the solution of which is 



where 



d_ 

dx 



dy(L/2) 
de 



2y(x) 



z(L/2) 



dy 

de 



z(x) = exp 



L/2 



" dx' 

L/2 Z(X'Y 



dx'y(x') 



Thus, Eq.(Al) can be written in the path-integral form 



(p(e)) = - \V'yG\y\e 



(A5) 



where S[y] is given by Eq.(A4), 



G[y] = ^{Lj2) 



+L ' 2 dx 

L/2 Z{x)' 



and the measure J D'y cx 53n=i ^Vn does not include the integrations over the values of y on the left and right 
boundaries. 

According to Eq.(|A3[), the optimal fluctuation 



U(x) 



2k 2 



cosh (kx) 

(see Eq.(p^|)) corresponds then to the instanton trajectory 

y(x) = Ktanh(Kcc), 



(A6) 



(A7) 



which describes the tunneling from the bottom of the l eft w ell, y — — n, to the bottom of the right well, y — +k and 
is the saddle-point trajectory of the Euclidean action (A4). Then the integration over disorder configurations close 
to the optimal fluctuation U reduces to the saddle-point integration of the imaginary-time trajectories close to the 
instanton trajectory y. The result of this integration is 

(p(e)) = ^rG[y]M-h-^ 1 , 



where 



M = det 



d 2 1 



In the limit large L, we obtain S[y] ~ |k 3 , G[y] ~ and the functional determinant M. can be found using the 
asymptotic behavior of the instanton solution Eq.(A7) at large distances£j 
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where the factor L comes from the integration over the instanton zero mode, which corresponds to the fact that the 
optimal fluctuation can be located at any place within the chain of size L. Combining these results, we obtain 

, , NX 4|e| 8H3/2 

{p{£)) = H e ~^ 



or, in usual units, 



which gives Eq. 



4lel 4v^R| £ | 3 / 2 

(pe = (A8) 



APPENDIX B: FIRST-ORDER CORRECTION TO THE WEIGHT OF THE OPTIMAL FLUCTUATION 

OF THE ELECTRON-HOLE PAIR 

The calculation of the first-order correction can be performed in two different ways. On the one hand we can use the 
exact relation between the action of the optimal fluctuation and the disorder potential averaged over the electron- hole 
wave function: 

d 2 x^ 2 ( Xl , x 2 ) {U{x x ) ~ 0U(x 2 )) = —-r fdxU 2 (x) = ~S. 



This equation allows us to write the action in the form 

A f, 9 , / h 2 d 2 h 2 d 



S=- d 2 X^ — - q — o + V(X! -X2)-E)V 



2 J \ 2m,h dx\ 2m e dx\ 

The first-order correction to the action is then given by 



S x = ^ Jd 2 x^> 2 V(x 1 - x 2 ) + SI, (Bl) 



where 

h 2 d 2 n 2 d 2 



S{ = A Jd 2 x^ 1 



2mh dx\ 2m e dx\ 



-E * 



'o 



Using Eq. (|34|) for the unperturbed electron-hole wave function <3>o and Eq.(|18|) for the single-particle wave function 
<fi, we obtain 

S( = 2\ d 2 x^ ^ 1 (\e h \(j) 2 (K h { Xl -x h )) + \e e \4> 2 (K e ( X2 - x e )j) . 



On the other hand, the direct calculation of the first-order correction to the action, in which one uses Eq.(^Tj), gives 
Sx = - £h f e [dx(/) 2 (K h (x - x h )<j) 2 (K e {x - x e )) + [dxU (x) Ui (x) , (B2) 





where U\ is the first-order correction to the disorder potential: 



U x {x) = -2A Jdx' [#q(x,x')$x(x,x') - /3^ q (x',x)^x(x',x)] 
The second term in Eq.(B2) can be rewriten in the form 

' dxU (x)Ux(x)=4:X Jd 2 x^o(xi,x2)^i(xi,X2) 

x \e h \4> 2 (K h (xx - Xh)) + \e e \(/) 2 (K e (x2 - x e )) (B3) 



-0 1 \£e\4> 2 (K e (x 1 ~ X e )) - j3\s h \(j) 2 (K h {x 2 - X h )) 



2S(, 
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where in the last step, we used that the overlap between the unperturbed electron and hole wave functions, iph{x — Xh) 
and ip e (x — x e ) is small, as a result the third and the fourth terms in the square brackets of Eq.( |B3|) give a much 
smaller contribution than the first and the second terms. Thus, we obtain 



Si 



4e/ 1 e e 



dxcj) 2 (Kh(x - x h ))<j) 2 {n e {x - x e )) + 2S{. 



Combining Eqs.(Bl) and (B4), we obtain 

Si = A Jd 2 x^lV( Xl 
which is equivalent to Eq.(p^). 



X2) 



^Shee 



dx(f) 2 {Kh(x - X h ))(f> 2 (He(x - X e )), 



(B4) 



(B5) 
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